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Q\ • Abstract 



The real, nonsingular elliptic solutions of the Korteweg-deVries equation are studied through 
the time dynamics of their poles in the complex plane. The dynamics of these poles is governed 
by a dynamical system with a constraint. This constraint is shown to be solvable for any 
finite number of poles located in the fundamental domain of the elliptic function, often in 
many different ways. Special consideration is given to those elliptic solutions that have a real 



^ , nonsingular soliton limit 



1 Introduction 

a' 



In 1974, Kruskal [15] considered the interaction of solitons governed by the Korteweg-deVries equa- 
tion (KdV), 

u t = 6uu x + u xxx . (1.1) 
Each KdV soliton is defined by a meromorphic function in the complex x-plane (i.e., sech 2 /c(x— xq)), 



so Kruskal [15] suggested that the interaction of two or more solitons could be understood in terms 
of the dynamics of the poles of these meromorphic functions in the complex x-plane, where the 
poles move according to a force law deduced from (|0|) . This was followed by the work of Thickstun 
|l7fl who considered the case of two solitons in great detail. 

Following a different line of thought, Airault, McKean and Moser Q studied rational and elliptic 
solutions of the KdV equation. An elliptic solution of KdV is by definition a solution of the KdV 
equation that is doubly periodic and meromorphic in the complex x-plane, for all time. Note that 
the soliton case is an intermediary case between the elliptic and the rational case. It was treated 
as such in 101. 
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Airault, McKean and Moser approached these elliptic solutions and their degenerate limits 
through the motion of their poles Xi(t) in the complex x-plane. In particular, they looked for elliptic 
KdV solutions of the form 

JV 

u(x,t) = -2j2p(x-Xi(t)). (1.2) 
i=i 

Here p(z) denotes the WeierstraB elliptic function. It can be defined by its meromorphic expansion 

P(Z) = I2+ £ {(z + 2mw x + 2nu 2 y ~ {2muj x + 2nu 2 )*) ' (L3) 

with u>x/u>2 not rea^]. More properties of the WeierstraB function will be given as they are needed. 
It is shown in || that the dynamics of the poles Xi(t) is governed by the dynamical system 

JV 

Xi = 12 p{xi~Xj), i = l,2,...,N, (1.4 a) 

(the dot denotes differentiation with respect to time) with the invariant constraint 

N 

p'(x l -x J ) = 0. i = l,2,...,N. (1.4 b) 



Here the prime denotes differentiation with respect to the argument. The solutions ( |l.2|) generalize 
an elliptic solution given earlier by Dubrovin and Novikov j7j, corresponding to the case N = 3. 
These authors also recall the Lame-Ince potentials [13] 

u(x) = -g{g + l)p(x), (1.5) 
which are the simplest <?-gap potentials of the stationary Schrodinger equation 

-^ + u{x)^ = \i,. (1.6) 

The remarkable connection between the KdV equation and the stationary Schrodinger equation 
has been known since the work of Gardner, Greene, Kruskal and Miura 11. Dubrovin and Novikov 
show @ that the (N = 3)-solution discussed in [fjj is a 2-gap solution of the KdV equation with a 
2-gap Lame-Ince potential as initial condition. 

If one considers the rational limit of the solution ( |1.2|) (i.e., the limit in which the Weier- 
straB function p(z) reduces to 1/z 2 ), then the constraint ( |1.4 b| ) is only solvable for a triangular 
number of poles, 

N = ^±», (1.7) 

for any positive integer n 0|. Notice that the Lame-Ince potentials are given by g(g J t- 1)/2 times an 
N = 1 potential. Based on these observations, it was conjectured in Q that also in the elliptic case 
given by (1.2), the constraint ( 1.4 b| ) is only solvable for a triangular number N, "or very nearly 

1 In this paper i t is always assumed that W\ is real and u>2 is imaginary. This is necessary to ensure reality of 
the KdV solution (1.2) when x is restricted to the real x-axis. Other considerations for reality of the elliptic KdV 
solutions will be discussed in Section [|. 
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so" . From the moment it appeared this conjecture was known not to hold, because it already fails 
in the soliton case, where the Weierstrafi function degenerates to hyperbolic functions. This failure 



of the conjecture easily follows from the work of Thickstun [17|. 

A further understanding of the elliptic case had to wait until 1988, when Verdier provided more 



explicit examples of elliptic potentials of the Schrodinger equation 21]. Subsequently, Treibich and 
Verdier demonstrated that 



u{x) 



4 

i=i 



9i{g% + 1) 



p(x -Xq- UJi) 



(1.8) 



the gi are positive integers) are finite-gap potentials of the stationary 



and hence result in elliptic solutions of the KdV equation IS, 18, 



(o>3 = OJl + UJ2,0Ji = 0, 
Schrodinger equation (JTT( 

The potentials of Treibich and Verdier were generalized by Gesztesy and Weikard [10, |llf . They 
showed that any elliptic finite-gap potential of the stationary Schrodinger equation ( |1.6| ) can be 
represented in the form 



M 



u(x) 



-2}^ 7, PK X 



a; 



(1.9) 



for some M and positive integers gi. Notice that this formula coincides with (|L2j ) if all the gi are 
1. 

The focus of this paper is the constrained dynamical system ( 1.4 a[ -b). We return to the ideas 
put forth by Kruskal [15] and Thickstun [[Pfl| . This allows us to derive the system ( |1.4 a| -b) in a 
context which is more general than |^|: there it was obtained as a system describing a class of 
special solutions of the KdV equation. Here, it is shown that any meromorphic solution of the KdV 
equation which is doubly periodic in x is of the for m (|1.2[ ). Hence the consideration of solutions of 
the form fll.2j ) and the system of equations ( 1.4 a - |l.4 b[) leads to all elliptic solutions of the KdV 
equation. 



Simultaneously, some of the results of Gesztesy and Weikard [1C, 11] are recovered here. Be- 
cause of the connection between the KdV equation and the Schrodinger equation, any potential of 
Gesztesy and Weikard can be used as an initial condition for the KdV equation, which determines 
any time dependence of the parameters of the elliptic KdV solution with that initial condition. Our 
approach demonstrates which parameters in the solutions Ql.9j ) are time independent and which 
are time dependent. 

The following conclusions are obtained in this paper: 



All finite-gap^] elliptic solutions of the KdV equation are of the form (|1.2| ) , for almost all times 
(see below). In other words, if u(x,t) is a finite-gap KdV solution that is doubly periodic in 
the complex x-plane, then u necessarily has the form ( |l.2j ) except at isolated instants of time. 



Any number N ^ 2 of Xj is allowed in ( |1.2|) , if the ratio |u;i/u;2| 3> 1. It follows that in 
this case, the constraint ( |1.4 bj ) is solvable for any positive integer^] N ^ 2. Our method 
for demonstrating this is to use a deformation from the soliton limit of the system ( |1.4 a[ -b). 
Furthermore, this method is useful because it provides good initial guesses for the numerical 
solution of (|1.4 b ) in Section ||. This method does not recover all elliptic solutions of the KdV 



2 See Section || 
3 The constraint 
for N = 1, with smaller periods 



1.4 b) is solvable for N = 2 [g]. As noted there, the corresponding solution reduces to a solution 
This case is therefore trivial and is disregarded 
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equations, because only those solutions are found which have nonsingular soliton limits. In 
particular, the solutions corresponding to the Treibich-Verdier potentials | ) are not found. 

If / <i^2 1 S> 1, then for a given N > 4, nonequivalent configurations satisfying the constraint 
exist that do not flow into each other under the KdV flow and which cannot be translated 
into each other. To the best of our knowledge this is a new result. 

The Xi are allowed to coincide, but only in triangular numbers: if some of the Xi coincide 
at a certain time t c , then gi{gi + l)/2 of them coincide at that time t c . At this time t c , the 
solution can be represented in the form ( |1.9[ ) with not all gi = 1. 

Such times t c are referred to as collision times and the poles are said to collide at the collision 
time. Before and after each collision time all Xi are distinct, hence pole collisions are isolated 
events. At the collision times, the dynamical system (|1.4 a|) is not valid. The dynamics of 
the poles at the collision times is easily determined directly from the KdV equation. 

Gesztesy and Weikard |l0|, |ll| demonstrate that ( |1.9| ) are elliptic potentials of the Schrodinger 
equation. These potentials generalize to solutions of the KdV equation, but this requires the 
gi to be nonsmooth functions of time. Only at the collision times t c are the gi not all one. 
Furthermore, at all times but the collision times, the number of parameters aj (which are 
time dependent) is N = X)i=i 9i (di + l)/2> This shows that the generalization from |ll|] 
to solutions of the KdV equation is nontrivial. 

Using the terminology of Chapter 7 of Q , the poles are in general position if all gi are equal 
to one. Otherwise, if not all gi = 1, the poles are said to be in special position. We conclude 
that at the collision times the poles are in special position. Otherwise they are in general 
position. 

The solutions discussed here are finite-gap potentials of the stationary Schrodinger equation 
(|l.6| ), with t treated as a parameter. In other words, each solution specifies a one-parameter 
family of finite-gap potentials of the stationary Schrodinger equation. It follows from our 
methods that to obtain a g-gap potential that corresponds to a nonsingular soliton potential, 
one needs at least N = g(g + l)/2 poles Xj. The Lame-Ince potentials show that this lower 
bound is sharp. If we consider potentials that do not have nonsingular soliton limits (such as 
the Treibich-Verdier potentials fll.SD ) then it may be possible to violate this lower bound. 



The relationship between elliptic solutions and soliton solutions of ( 1.4 a| -b) and hence of the 



KdV equation is made explicit. In fact, this relationship is essential to the method used here, 
as mentioned above. 

In Section (||), we present an explicit solution of the form ( |l.2j ) with N = 4. Notice that if in 
Ql.SD all gi are one, the solution is reducible to a solution with N = 1 and smaller periods. To 
see this it is convenient to draw the pole configuration corresponding to ( |l.8| ) in the complex 



plane. This is actually true, even if gi is not one, but all gi are equal. In that case, (1. 
reduces to a Lame-Ince potential ( |1.5[ ). Unlike any of the Treibich-Verdier solutions, the 
(N = 4)-solution, presented in Section (||), has a nonsingular soliton limit. 

The first five conclusions are all discussed in Section In Section |2j, the results of Kruskal 
|l5|| and Thickstun (lTj for the soliton solutions of the KdV equation are reviewed, but they are 
obtained from a point of view which is closer to the approach we present in Section || for the elliptic 
solutions of the KdV equation. Finally, in Section [|, some explicit examples are given, including 
illustrations of the motion of the poles in the complex x-plane. 
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2 The soliton case: hyperbolic functions 



In this section, the results of Kruskal |15| and Thickstun [17] for the dynamics of poles of soliton 
solutions are discussed from a point of view that will allow us to generalize to the periodic case. 
Consider the one-soliton solution of the KdV equation 

u(x, t) = 2k 2 sech 2 k(x + Ak 2 t - <p). (2.1) 

Here A; is a positive parameter (the wave number of the soliton) determining the speed and the 
amplitude of the one-soliton solution. Using the meromorphic expansion [|E|] 

1 x 00 1 

—rrcosech 2 — = —77, (2.2) 

(uniformly convergent except at the points x = innT) one easily obtains the following meromorphic 
expansion for the one-soliton solution of the KdV equation: 



U(x,t) = -2k 2 V — — 5 r 777- (2.3) 

v ' ; n ^ jo (k(x + Ak 2 t - <p) + if + mvr) 2 v ; 

From this expression, one easily finds that the locations of the poles of the one-soliton solution of 
the KdV equation for all time are given by 

2 , VK ( 1 



x n = if - 4k z t - — I n + - ) . ( 2 ' 4 ) 



This motion is illustrated in Fig. 2.1a. Notice that the locations of the poles are symmetric with 
respect to the real x-axis. This is a consequence of the reality of the solution (|2.l| ). In order for 
a solution to be real it is necessary and sufficient that if x n (t) is a pole, then so is x* n (i), where 
* denotes the complex conjugate^. The closest distance between any two poles is d = ir/k and is 
constant both along the vertical line Re(x) = (p — 4k 2 t and in time. Note that the poles are moving 
to the left. This is a consequence of the form of the KdV equation (|1.1|), which has time reversed, 
compared to the version Kruskal [15] and Thickstun [17 used. 

Since two solitons of the KdV equation cannot move with the same speed, a two-soliton solution 
of the KdV equation asymtotically appears as the sum of two one-soliton solutions which are well- 
separated: the higher-amplitude soliton, which is faster, is to the right of the smaller-amplitude 
soliton as t — > — oo. Ast^oo, the higher- amplitude soliton is to the left of the smaller-amplitude 
soliton. Hence as t — * — oo, the pole configuration of a two-soliton solution with wave numbers k\ 



and &2 is as in Fig. |2.1[ b. In this limit, the two-soliton solution is a sum of two one-soliton solutions. 
Each results in a vertical line of equispaced poles, with interpolar distance respectively d\ = n/ki 
and d2 = vr//c2- As long as the solitons are well-separated, these poles move in approximately 
straight lines, parallel to the real axis, with respective velocities v\ = —Ak 2 and V2 = — 4fc|- Since 
\ v i\ > | ^2 1 5 the solitons interact eventually. This interacting results in non-straight line motion of 
the poles. After the interaction, the situation is as in Fig. [2.1| b, but with the two lines of poles 
interchanged. 



Thickstun [17| considered the case where k\ and &2 are rationally related, so kxjki = p/q, where 



p and q are positive integers. In this case, one can define D = pd\ = qd2- The complex x-plane 



4 The vertical line of poles can be rotated arbitrarily. The expression (2.3) still results in a solution of the KdV 
equation, but it is no longer real. Again, we will only consider real, nonsingular solutions, when restricted to the real 
a;- axis 
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v=-4 k 



D 



v,=-4 k , 



1 



d= 



v,=-4 1 



(a) 



(b) 



t 



Figure 2.1: (a) The motion of the poles of a one-soliton solution in the complex plane, 
(b) The asymptotic motion of the poles of a two-soliton solution. 



is now divided into an infinite number of equal strips, parallel to the real x-axis, each of height 
D. The real x-axis is usually taken to be the base of such a strip. It is easy to show )i~7| ] that the 
motion of the poles in one strip is repeated in every strip. Hence, one is left studying the motion 
of a finite number N (= p + q) of poles in the fundamental strip, whose base is the real x-axis. 
Thickstun examined this motion by analyzing the exact expression for a two-soliton solution of the 
KdV equation. 

Any two-soliton solution is expressible as [Q] 

u(x,t) = 2dllnT{x,t). (2.5) 

It follows from this formula that the poles of u(x,t) are the zeros of r(x,t) if r(x, t) is entire in x. 
Then the Weierstrafi Factorization Theorem || gives a factorization for r(x,t): 



T(X, 

k=l 



'*) =C Ii I 1 )e x ' Xk . (2.6) 



Since only the second logarithmic derivative of this function is relevant, the constant C is not 
important. If the solution is periodic in the imaginary x-direction, this is rewritten as 



N oo 



rOM) = Cn II ( 1 ~ ~^~T7rn ) e x ' ( - Xn+ilD \ (2.7) 

n=l l=— oo 



where the first product runs over the poles in the fundamental strip. The second product runs over 
all strips. Using the uniform convergence of (12.71) > 

N oo 

" (MH - 2 S,£(^w (2 - 8> 
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which, using (2.2), is rewritten as 



2 N 



u(x,t) = —2—— cosech 



D 



k 2 p q 
-2-| ^2 cosech z 



k 2 (x 



(2.9) 



=1 " ~* n=l 

where the pole locations x n depend on time: x n = x n (t). One recovers the one-soliton solution ( |2,lD 
easily, by equating k\ = 0,p = 0, q = 1. Equation ( [2.9P essentially expresses a two-soliton solution 
as a linear superposition of iV one-soliton solutions with nonlinearly interacting phases. Note that 
the first equality in ( |2.9| ) is valid for arbitrary soliton solutions that are periodic in x with period iD. 
This is the case for a g-soliton solution if its wavenumbers hi, i = 1, 2, . . . , g are all commensurable: 
(ki : k 2 : . . . : k g ) = (pi : p 2 : . . . : p g ), for positive distinct integers pi, i = 1, 2, . . . , g which have no 
overall common integer factor. The total number of poles in a strip is then N = p\ +P2 + • • • + Pg- 
In obtaining (E 



and fl2.S|), we have deviated from Thickstun's approach [17| to an approach that 
is generalized to the elliptic case of the next section in a straightforward way. 

Next, we derive the dynamics imposed on the poles x n (t) by the KdV equation. This is con- 
veniently done by substituting ( |2.£| ) into fll.ip and examining the behaviour near one of the poles: 
x = x n + e. This results in several singular terms as e — ► 0, corresponding to negative powers of e. 
The dynamics of the poles is then determined by the vanishing of the coefficients of these negative 
powers and the zeroth power. This results in only two nontrivial equations, obtained at order e~ 3 
and e -2 respectively: 



for n = 1, 2, 







N 



12 E E 

N oo 

E E 



i 



(x k - x n - iwy 



+ 12 J2 



i^0,/=-oo 



1 

(x k - x n - UDy 



k=l,k^n l=—oo 

, N. Using (|2.2|) and its derivative, 



(2.10 a) 
(2.10 b) 



N 



-4— j + 12— j > coseclr 
D 2 D 2 t-i, 

k=l,k^en 



ir(x k 



N 



cosech 2 

k=l,k^n 



7T(x k - X n ) 

D 



coth 



D 

n(x k - x r 



D 



(2.11 a) 



(2.11 b) 



for n = 1, 2, . . . , N. Hence the dynamics of the poles x n (t) is determined by ( 2.11 a ). This dynamics 
is constrained by the equations ( [2.11 b ). These constraint equations ( 2.11 b| ) are invariant under 
the flow of ( 2.11 a| ). This follows from a direct calculation. 



Remarks: 

• Since the KdV equation has two-soliton solutions for any ratio of the wavenumbers k\/k 2 ^ 1, 
the constraint ( |2.11 b ) is solvable for any value of N, excluding N = 2, which can only be 
obtained by p = q = 1, resulting in equal wavenumbers k\ and k 2 - 

• In particular, it follows that the minimum number of poles in a fundamental strip required 
to obtain a g — soliton solution isiV = 1 + 2 + 3 + .. . + g = g(g + l)/2, corresponding to 



a g-soliton solution with wavenumbers which are related as 
.'/ 1:...:2:1). 



(h : k 2 



'■ k n 



-l 



(9-- 
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• Equating k\ = 0, p = 0, q = 1, one obtains from ( 2.11 a| ) x\ = —4^1, corresponding to the 
dynamics of the one-soliton solution. The asymptotic behavior of the poles of a two-soliton 
solution also follows from ( |2.11 a ): from the separation of the poles into distinct vertical 
lines, it follows from ( 1.4 a| ) that the velocity of these vertical lines is given by the one-soliton 
velocity for each line, as expected. This result follows from easy algebraic manipulation and 
the identity 

n 2 1 P- 1 



P 



1 



E 

n=l 



cosec 



nir 



P 



(2.12) 



valid for any integer p > 1. 



A full analysis of the interaction of the poles for the case of any two-soliton solution with 
=p/q is given in [17]. 



3 The elliptic case 



Consider the quasiperiodic finite-gap solutions of the KdV equation with g phases [14| 

u(x,t) = 2dlln6 g (kx + ut + <j>\B), 



where 



6 g (z\B) 



£ exp 
meZ 9 



1 



-m ■ B ■ m + im ■ z 



(3.1) 



(3.2) 



a hyperelliptic Riemann-theta function of genus g. The g x g real Riemann matrix (i.e., symmetric 
and negative definite) B originates from a hyperelliptic Riemann surface with only one point at 
infinity. Furthermore, k, u and (j) are (7-dimensional vectors. 

The derivation of equations ( |2.9| ), (2.11 a) and ( [2.11 b| ) is easily generalized to the case where 
the solution is not only periodic in the imaginary x-direction, but also in the real x-direction: 



u(x + Li,t) = u(x, t) = u(x + iL,2,t). 



(3.3) 



This divides the complex x-plane into an array of rectangular domains, each of size L\ x L2. One 
of these domains, called the fundamental domain S, is conveniently placed in the lower left corner 
of the first quadrant of the x-plane. The theta function has the property H 



g (z + iB ■ M + 2ttN\B) = 9 g (z\B) exp 



-M ■ B ■ M + iM ■ z 

2 



(3.4) 



for any pair of g-component integer vectors M, N. This expression is useful to determine conditions 
on the wavevector k in order for u(x,t), given by (3.1), to satisfy (3.3): 



3 No, M Q £ Z 9 : kL x = 2irN , kL 2 = B ■ M . (3.5) 

These results are now used to determine the number of poles N of u(x, t) in the fundamental 
domain. The poles of (3.1) are given by the zeros of i?(x,t) = 9 g (kx + ut + 4>\B), regarded as a 
function of x: 
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N = — I d\n$(x,t) 

JOS 



1 / fL! pL 2 

/ dh\-d(x,t) + / dlntf(x + L x ,t)- 

2m yjo Jo 

J dlntf(x + iL 2 ,t) - J dln$(x,t)\ 



using (3.4) 1 f L i , ti(x,t) 
— ' a In 



2m' Jo + iL2,t) 

using (p^i ) 1 Z" 1 ! 



2vri 



/ d(-iiM • fe) 
Jo 



-M -N = --^-M -B-M . (3.6) 

The first equality of ( |3.6| ) confirms that A?" is an integer. The second equality shows that N is 
positive, by the negative-definiteness of B. 

We now proceed to determine the dynamical system satisfied by the motion of the iV poles of 
u(x,t) in the fundamental domain S. Again, the poles of u(x,t) are the zeros of i?(x,t). Further- 
more, simple zeros of $(x, t) result in double poles of u(x, t), as in the hyperbolic case. 

The Weierstrafi Factorization theorem || gives the following form for i}(x,t): 

&(x, t) = e cx2 ' 2 nfl--) e^ + ^l , (3.7) 



where the product runs over all poles xt- The additional exponential factors, as compared to 
( |2.6| ), are required because the poles now appear in a bi-infinite sequence: both in the vertical and 
horizontal directions. These exponential factors ensure uniform convergence of the product. The 
parameter c is allowed to depend on time. It determines the behavior of $(x, t) as x approaches 
infinity in the complex x-plane ||. Using (|3.3|), this is rewritten as 



N oo oo 

^(m) = ex P (cx 2 /2) n n n i 1 

n=l m=— oo l=~oo 



x n + mil + UL 2 



o 

X X 



6XP y x n + mil + UL 2 + 2{x n + mLi + UL 2 ) 2 J ' 

The first product runs over the number of poles (N) in the fundamental domain, the second and 
third products result in all translations of the fundamental domain. From the uniform convergence 
of (O), 



N oo oo r 

(x,t)=2c-2j2 E E ( (x _ Xn _ mLl _ ilL2) 2 - {Xn + mLl + llh2) 2 ) ■ ( 3 - 9 ) 



n=l m=—oo l=—oo 



Using the definition of the Weierstrafi function ( |1.3[ ), this is rewritten as 
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N N 

u(x, t) = 2c - 2^2 p(x - x n ) + 2^2 p{x n ), 
n=l j=l 

where the periods of the Weierstrafi function are given by 2u\ = L\, 2oj2 = ^2- Define 



(3.10) 



N 



c = 2c + 2^p(x n ). 



(3.11) 



n=l 



The dynamics of the poles x n = x n (t) is determined by substitution of ( [3.10 ) or ( |3.9[ ) into the KdV 
equation and expanding in powers of e for x near a pole: x = x^ + e. Equating the coefficients of 
e -3 , e -2 and e° to zero result in 



N 

12 £ p( Xj 

AT 



j=l,j^n 

c(t) = a = 0, 



(3.12 a) 

(3.12 b) 
(3.12 c) 



for n = 1,2, . . . , N . (The constant a can always be removed by a Galilean shift, so it is equated 
to zero, without loss of generality.) The constraints (|3.12 b| ) are invariant under the flow, as can 
be checked by direct calculation. Notice that ( 3.12 a[ -b) are identical to the equations obtained 
by Airault, McKean and Moser Q. These equations are obtained here in greater generality: any 
solution Q3.1| ) which is doubly periodic in the x-plane gives rise to a system ( 3.12 a| -b). This allows 
us to reach the conclusions stated in the next section. 



Remarks 



In the limit L\ — > oo, the equations (3.12 a- 3.12 b| ) reduce to ( [2.11 a -2.11 b). This limit is 
most conveniently obtained from the Poisson representation of the Weierstrafi function: 



it \ 2 1 2 irx ^ r .o 71 "/ t \ , 2 n7rLi 

— I - + cosecn — — h > < cosecn — (x + nL\) — cosecn — - — 
L2J \ 3 L2 , n L L2 L2 



(3.13) 

This representation is obtained from ( |1.3[ ) by working out the summation in the vertical 
direction. It gives the Weierstrafi function as a sum of exponentially localized terms, hence 
few terms have important contributions in the fundamental domain. A Poisson expansion for 
p'{x) is obtained from differentiating ( p.!3| ) term by term with respect to x. 



• Define the one-phase theta function 9i{z,q) |12j : 



9i(z,q) = 2^(-l) r y n+1 / 2 ) 2 sin(2n+ l)z, 

n=0 



(3.14) 



with \q\ < 1. If L2 < L\, then the relationship p(z) = a — d^ln6i(Trz/Li,iL2/Li), with a a 



constant [12], allows us to rewrite (1.2) as 
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N , r \ 

u(x,t) = a + 2d 2 x In TT 6>i ( -^(x - &,-(*)), i-^ ) , (3.15) 

with a = —2aN. Hence, for the doubly-periodic solutions of the KdV equation of the form 
Q3.1| ), it is possible to rewrite the g-phase theta function as a product of N 1-phase theta 
functions, with nonlinearly interacting phases. Note that this does not imply that the g- 
phase theta function appearing in ( |3.1| ) is reducible. Reducible theta-functions do not give 
rise to solutions of the KdV equation H. 

By taking another time derivative of ( |3.12 a| ) and using ( |3.12 b| ), one obtains 

N 

x n = -(12) 2 P'( X j ~ x n)p{Xj - x n ). (3.16) 

It is known that this system of differential equations is Hamiltonian Q , with Hamiltonian 

, N / 19 n2 N N 

H = \Y,Pl+ { -^-Y, E f?i*k-xj), (3-17) 
k=i k=ij=i,j^k 

and canonical variables {xk,Pk = &k}- A second Hamiltonian structure for the equations 
( 3.16 ) is given in 0. A Lax representation for the system ( 3.12 a| -c) is also given there. This 



Lax representation is a direct consequence of the law of addition of the Weierstrafi function 
[12]. It is unknown to us whether a Hamiltonian structure exists for the constrained first-order 
dynamical system ( |3. 12 af b). 

The Hamiltonian structure (3.17) shows that the system ( [3,12 af b) is a (constrained) member 



of the elliptic Calogero-Moser hierarchy 



4 Discussion of the dynamics 



In this section, the constrained dynamical system ( |3.12 a -c) is discussed. In particular, the asser- 
tions made in the introduction are validated here. 



For reality of the KdV solution ( |3.lD when x is restricted to the real line, it is necessary and 
sufficient that if Xj(t) appears, then so does Xj(t). Because the Weierstrafi function is a meromorphic 
functions of its argument, this reality constraint is invariant under the dynamics ( 3.12 a| ). 



As a consequence, the distribution of the poles in the fundamental domain S is symmetric with 
respect to the horizontal centerline of S. Poles are allowed on the centerline. Most of what follows 
is valid for both real KdV solutions^] and KdV solutions that are not real, but we restrict our 
attention to real KdV solutions. 



4.1 All finite-gap elliptic solutions of the KdV equation are of the form ( |1.2| ), 
up to a constant 

A straightforward singularity analysis of the KdV equation 0] shows that any algebraic singularity 
of a solution of the KdV equation is of the type u(x,t) = —2/(x — a(t)) 2 + 0(x — a), for almost 

5 "Real KdV solution" refers to a solution of the KdV equation which is real when x is restricted to the real a;-axis 
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all times t. At isolated times t c , the leading order coefficient is not necessarily —2. It can be of the 
form —g(g + 1) (see below), but the exponent of the leading term is always —2. 

Hence, an elliptic function ansatz for u(x, t) can only have second order poles and with the 
substitution u{x,t) = 2d x \wd(x, t) gives rise to a Weierstrafi expansion of the form (|3.8|), with an 
arbitrary prefactor exp(c(x, t)), for an arbitrary function c(x,t), entire in x. Substitution of this 
ansatz in the KdV equation then determines that c xx is doubly periodic and meromorphic in x. 
The only such c is a constant. Hence, all finite-gap elliptic solutions of the KdV equation are of 



the form (1.2). 



4.2 If Li/L 2 ^> 1, any number of poles N ^ 2 in the fundamental domain is 
allowed 



We have already argued that the equations ( |2.9|) an d ( 2.11 a) -b) are obtained from ( |1.2p and ( p. 12 a 



b) in the limit L\ — > oo. On the other hand, ( |3.13|) can be regarded as a perturbative expansion of 
the soliton case, which corresponds to its first two terms. Because ( |3.13D converges uniformly away 
from x = n\Li + in 2 L2, for arbitrary integers rti,n2 we conclude that for large (but finite) values 
of L1/L2 any real, nonsingular (/-soliton solution with rationally related wave numbers (k% : k 2 : 
. . . : kg) = (pi : P2 ■ ■ ■ ■ '■ p g ) has an elliptic deformation with real period L\ and imaginary period 
1L2 , which is real and nonsingular^]. Since ./V = p\ + p 2 + . . . + p g 7^ 2 is arbitrary in the soliton 
case, this is also true for these elliptic deformations of the solitons. Hence the constraint ( [3.12 b ) 
is solvable for arbitrary N, for L\ 3> 

At this point, it is appropriate to remark that if one is interested in elliptic solutions of the 
Kadomtsev-Petviashvili (KP) equation, 

d x {~u t + 6uu x + u xxx ) + 3a 2 u yy = (4.1) 
(with parameter a), then the equations fll.2j ), ( |3. 12 a[ -b) are replaced by § 

N 

u{x,y,t) = -2^2p(x - x n (y,t)), (4.2 a) 

n=l 

2 N 



dXn 

dt 



S8r \ 2 

* a \dl +12 S P(^-^)' ( 42b ) 



d 2 x 



2 u *>n 

a 



dy 2 



16 P'{x 3 -x n ). (4.2 c) 



This clarifies the appearance of the constraint ( 3.12 b| ) on the motion of the poles of elliptic solutions 



of the KdV equation, where the poles are independent of y. The search for y-independent solutions 



of the KP equation reduces it to the KdV equation and it reduces the equations ( 4.2 a -c) to ( |l.2| 



( 3.12 a[ -b), forcing the poles to remain on the invariant manifold defined by ( 3.12 b| ). For the KP 



equation, no such constraint exists and hence the number of poles in the fundamental domain 
can be any integer, not equal to two. 



6 That this deformation is nonsingular follows from the fact that the limit of such a deformation is the original 
soliton solution: the limit of a singular solution results in a singular soliton solution. This is impossible, because 
only nonsingular soliton solutions are considered, hence the elliptic deformations of nonsingular soliton solutions are 
nonsingular 
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4.3 For any iV > 4, nonequivalent configurations exist, for L\/L 2 sufficiently 
large 

Again, we only consider solutions with nonsingular soliton limits; i.e., the elliptic deformations 
mentioned above. Consider the asymptotic behavior for t — ► — oo of lim^ 1 ^ 00 u(x, t). In this 
soliton limit, as t — > — oo, the poles are collected in groups corresponding to one-soliton solutions. 

In this section, two configurations are called nonequivalent if the above asymptotic behavior 
results in a different grouping of the poles. 

For N = 3, all configurations are equivalent to one configuration. In the limit L\ — > oo, this 
configuration corresponds to the two-soliton case with k\ : k 2 = 2 : 1. This configuration is 



discussed in Section 5.1 



For TV = 4, all configurations are again equivalent to one configuration. This configuration 
corresponds to the two-soliton case with k\ : k 2 = 3 : 1. Recall that k\ and k 2 are not allowed 
to be equal, hence a configuration with two poles to the left and two poles to the right does not 
exist. Another way of expressing that only one configuration exists is that TV = 4 can only be 
decomposed as the sum of distinct positive integers without common factor as N = 3 + 1. Again, 



all TV = 4 configurations are equivalent. This configuration is discussed in Section 5.2. That section 



also discusses another example of an N = 4 potential which does not have a nonsingular soliton 



limit. This potential is a special case of one of the Treibich-Verdier potentials (1.8). 

Any integer N > 4 can be written as a sum of distinct positive integers without overall common 
factor in more than one wayQ. Let the number of terms in the m-th decomposition of N be 
denoted as N m , then N = J2iJ^l n i> with the distinct and having no overall common factor. 
This configuration corresponds to the iV m -soliton case with wavenumber ratios {k\ : k 2 : . . . : 
kpf m ) = (ni ■ n,2 riN m ). A solution with these wave numbers has N m phases and is an 

iV m -soliton solution. Hence for any N > 4 there exist at least as many different configurations as 
there are decompositions of ./V into distinct positive integers, without overall common factor. These 
configurations need not have the same number of phases. 



Two nonequivalent configurations corresponding to N = 5 are discussed in Section 5.3 



4.4 The poles only collide in triangular numbers 

A collision of poles is a local process in which only the colliding poles play a significant role. The 
analysis of the collisions is identical to the rational and the soliton cases because close to the 
collision point, the WeierstraB function reduces to 1/x 2 . Kruskal |l5|] already noticed that the poles 
do not collide in pairs, but triple collisions do occur. In fact, any triangular number of poles can 
participate in a collision, in which case the solution of the KdV equation at the collision time t c , 
nearby the collision point x c is given by u(x,t c ) = —g{g + l)/(z — x c ) 2 Asymtotically near the 
collision point x c , i.e., t < t c , the poles lie on the vertices of a regular polygon with g(g + l)/2 
vertices. For t > t c , the poles emanate from the collision points, again forming a regular polygon 
with g(g + l)/2 vertices. If g(g + l)/2 is even, this polygon is identical to the polygon before the 
collision. If g(g + l)/2 is odd, the polygon is rotated around the collision point by 2n/g(g + 1) 
radians. 

Of all these collision types, the one where three poles collide (corresponding to g = 1) is generic. 
It is the one observed in the examples illustrated in Section |5[ 

Since the poles only collide in triangular numbers, it is possible that at any given time t c the 
solution of the KdV equation has the form ( |1.9[) , with not all g\ = 1. At almost every other time 

7 5 = 1 + 4 = 2 + 3, N = 1 + (N - 1) = 1 + 2 + (TV - 3), for TV > 5 
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t, such a solution has N = J2f=i 9i{di + l)/2 distinct poles. 

4.5 The solutions ( |1.2|) are finite-gap potentials of the stationary Schrodinger 
equation ( |1.6|) 



By construction the solutions ( |1.2| ) are periodic in x because they are obtained as a Weierstrafi fac- 
torization of the theta function appearing in (|3,1|), upon which we have imposed the double pe- 



riodicity. Hence the solutions (1.2) are finite-gap potentials of the Schrodinger operator. In |iT| , 
another proof of this can be found. 

For solutions that are elliptic deformations of the nonsingular solitons of Section |2], more can 
be said: an elliptic deformation of a g-soliton solution is a g-gap potential of the Schrodinger 
equation. The reasoning is as follows: we already know that any elliptic deformation results in a 
finite-gap potential of the Schrodinger equation. On the other hand, any finite gap potential of 
the Schrodinger equation is of the form (|3.1|). The soliton limit of such a finite-gap solution with 
g-phases is a g-soliton solution ||. Hence the number of phases of an elliptic deformation of a 
g-soliton solution is equal to g. 

This limit is the soliton limit of the periodic solutions, in which the fundamental domain reduces 
to the fundamental strip. In order to have a g-soliton solution of the KdV equation, we remarked 
in Section |2| that at least N = g(g + l)/2 = 1 + 2 + .. . + (7 poles are required in the fundamental 
strip. Hence, this many poles are required in the fundamental domain to obtain a <?-gap potential 
of the Schrodinger equation that is an elliptic deformation of a nonsingular soliton solution. 



5 Examples 

In this section, some explicit examples of elliptic solutions of the KdV equation are discussed. These 
are illustrated with figures displaying the motion of the poles in the fundamental domain. Other 
figures display the solution of the KdV equation u(x, t) as a function of x and t. All these figures 
were obtained from numerical solutions of the corresponding constrained dynamical system. In 
all cases, the constrained dynamical system was solved using a projection method: the dynamical 
system ( |3.12 a ) is used to evolve the system for some time. Subsequently, the new solution is 



projected onto the constraints ( [3.12 b ) to correct numerical errors, after which the process repeats. 



In all examples given, L\ = 4 and L2 = vr. This seems to indicate that one can wander far 
away from L1/L2 3> 1 and still obtain soliton-like elliptic solutions of the KdV equation. This is 
not surprising as ( |3.13| ) indicates that as perturbation parameter on the soliton case one should use 
e = exp(— 2ttL\/L2). For the values given above, this gives e = 0.00034. 

5.1 The solution of Dubrovin and Novikov N = 3 

Dubrovin and Novikov integrated the KdV equation with the Lame-Ince potential u(x, 0) = 
— 6p(x — x c ) as initial condition. They found the solution to be elliptic for all time, with N = 3. 
They gave explicit formulae for the solution, which they remarked was probably the simplest 
two-gap solution of the KdV equation. The dynamics of the poles in the fundamental domain 



is displayed in Fig. 5.1a. Fig. 5.1b displays the corresponding two-phase solution of the KdV 
equation. Animations of the behavior of the poles and of u(x, t) as t changes are also available at 
trttp : //amath-www . Colorado . edu/appm/other/kp/papers| . Notice the soliton-like interactions of 



the two phases in the solution. In terms of the classification of Lax [16], these are interactions of 
type (c) (i.e., u(x,t) has only one maximum while the larger wave overtakes the smaller wave). 
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(a) (b) 

Figure 5.1: The solution of Dubrovin and Novikov, with L\ = 4 and L2 = vr. (a) The 
motion of the poles in the fundamental domain. The initial position of the poles is 
indicated by the black dots. The arrows denote the motion of the poles, (b) The KdV 
solution u(x, t) 



From Fig. 5.1a, it appears that the Dubrovin-Novikov solution is periodic in time. This was 
indeed proven by Enol'skh ||]. 

For this specific solution only one of the three constraint equations is independent: since the 
derivative of the Weierstrafi function is odd, the sum of the constraints is zero. Furthermore, 
labelling the three poles by xi,X2 and X3, for reality X2 = x\ + 1L2 and £3 is on the centerline. 
Hence the second constraint is the complex conjugate of the first constraint. The constraints 
( 3.12 b| ) reduce to the single equation 



p'(x 1 - x\) + p'(xi - x 3 ) = 0. (5.1) 



This equation was solved numerically to provide the initial condition shown in Fig. p.l|a. The 
initial guess required for the application of Newton's method is based on the knowledge of the 
soliton limit. In that case two poles on the right represent a faster soliton, one pole on the left 
represents the slower soliton. The periodic case is not that different: the vertical line of poles with 
the smallest vertical distance between poles has poles closer to the real x-axis than the others and 
correspond to the wave crest with the highest amplitude, as seen in Fig. [5.1|b. We refer to the 
Dubrovin-Novikov solution as a (2 : l)-solution because of the natural separation of the poles in a 
group of 2 poles (xi,X2) and a single pole (X3). 

Equating x\ = x$ + e and only condering the singular terms of (|5.1|), it is possible to examine the 
location of the poles close to a collision points x c . With p'(x) = —2/x s in this limit and e = e r + iei, 
one finds 

el - 3ef e r = 0, 9ef - 3e^ = 0. (5.2) 
This set of equations has three solutions, corresponding to the three distances between the poles: 



15 



e r /ei G |o, \/3, — \/3|. This allows for two triangular configurations of the poles: an equilateral 
triangle pointing left of the collision point and one pointing right. 

Using the dynamical system ( 0.12 a| ) in the same way and only retaining singular terms results 



m 



■T2+ 3 - 



6e r e, ; 



(5.3) 



Since the constraints ( 0.12 b ) are invariant under the flow (|3.12 a[ ), the solutions to ( |5.2[ ) give 
invariant directions of the system ( |5.3j ). Along these invariant directions, one obtains ordinary 
differential equations for the motion of the poles as they approach the collision point. It follows 
from these equations that the poles approach the collision point x c with infinite velocity. Integrating 
the equations with initial condition e(t c ) = gives 



3 5/6 



(V3 + i)(t c -t 



,V3 



(5.4) 



Using the three branches of (t c — i) 1 ' 3 results in the dynamics of each edge of the equilateral triangle. 
If t < t c this triangle is pointing left, for t > t c it is pointing right. 

5.2 N = 4: an elliptic deformation and a Treibich-Verdier solution 

The next solution we discuss has 4 poles in the fundamental domain and is an elliptic deformation 
of a soliton solution. In the limit L\ — * oo, this solution corresponds to a two-soliton solution with 
wavenumber ratio fci/fo = 3/1, so this solution is refered to as a (3 : 1)- solution. 



3 
2.5 

2 
1.5 

1 

0.5 



0.5 1 1.5 



2.5 



3.5 



Figure 5.2: = 4: the pole dynamics of a (3 : l)-solution. The initial positions of the 4 
poles are indicated. The arrows on the centerline indicate that the poles there move 
in both directions. 



The motion of the poles in the fundamental domain is displayed in Fig. 5.2. Corresponding 



to the given wavenumber ratio, the amplitute ratio of the two phases present in the solution is 
roughly k\jk\ = 9/1. As a consequence, the form of u(x, t) is not very illuminating and it has been 
omitted. Animations with the time dependence of both the positions of the poles and of u(x, t) are 
again available at http : //amath-www. Colorado . edu/appm/other/kp/papers| . 
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Note that the poles of the (3 : l)-solution do not collide. This is in agreement with the results 



of Thickstun [17] who outlined which configurations lead to collisions and which do not, in the 
hyperbolic case. As mentioned before, the examination of collision behavior is essentially local and 
no differences appear between the rational, hyperbolic and elliptic cases. 
Another configuration with N = 4 exists. Consider the potential 

u(x, t = t c ) = —2p(x — xq) — 6p(x — xq — &x), (5-5) 



with xq on the centerline. This is a Treibich-Verdier potential, obtained from (1.9) with M = 2, 
9i = 1) 9i = 2, a\ = xo and a-z = Xq + u%. It is referred to as a Treibich-Verdier potential because 
the position of the poles is given in terms of the periods of the Weierstrafi function, as in (|1.8| ). Also, 
it can be obtained from fll.S| ) as a degenerate case. As before ./V = Yli= l 9i(di + l)/2 = 4, hence for 
all times that are not collision times, this solution has 4 distinct poles in the fundamental domain. 
The time t = t c is a collision time. Immediately after the collision time t = t c , the 3 poles located 
at xo + oj\ separate, as in the Dubrovin-Novikov solution, along an equilateral triangle. The result 
appears to be a three-phase solution. However, it is known that the potential (|5.5D is a two-gap 
potential of the Schrodinger equation and its hyperelliptic Riemann surface is given explicitly in 
H . This solution is not an elliptic deformation of a nonsingular soliton solution and the separation 
into different phases does not make sense. This is also seen from the following argument: if, for a 
fixed time which is not a collision time, we attempt to take the limit as L\ = 2ui — ► oo, the poles 
seem to separate in three distinct solitons with respective wave numbers (ki : k% '■ k%) = (2:1:1). 
Such a nonsingular soliton solution does not exist for the KdV equation and the separation into 
different phases does not make sense. 



The dynamics of the poles is illustrated in Fig. |5.3| a. The corresponding KdV solution is shown 
in Fig. |Jb. 

The dynamics of the poles illustrated in Fig. |5li| a exhibits behavior that appears qualitatively 
different from any other solution discussed here. The trajectories traced out by the motion of the 
poles in the fundamental domain appear to have singular points (cusps), away from the collision 
points. Upon closer investigation, these "cusps" are only a figment of the resolution of the plot and 
the poles trace out a regular curve as a function of time, away from the collision times. Exactly 
why the global pole dynamics of the Treibich-Verdier potential (|5.5|) under the KdV flow appears so 
different from the pole dynamics of elliptic deformations of soliton solutions of the KdV equation is 
an open problem. Another question one may ask is whether similar behavior is observed for other 
solutions originating from Treibich-Verdier potentials. 

5.3 N = 5: two different possibilities 

For N = 5, two soliton configurations are possible, and corresponding to each of these is an elliptic 
solution. The first solution is a (4 : l)-solution. The second solution is a (3 : 2)-solution. 

The (4 : l)-solution offers no new pole-dynamics: initially 1 pole is located on the centerline, 
at the left in the fundamental domain. The other poles are located at the right of the fundamental 
domain, symmetric with respect to the centerline. The three poles closest to the centerline interact 
as the (2 : l)-solution. The two outer poles behave as the two outer poles of the (3 : l)-solution. 



The pole dynamics of the (4 : l)-solution is displayed in Fig. 5.4a. 



The (3 : 2)-solution is more interesting. It is displayed in Fig. |5.4|c, together with the motion of 



the poles in the fundamental domain 5.4b. Again, the two crests of u(x,t) interact in a soliton-like 



manner. In Lax's classification [16], this is an interaction of type (a), where at every time two 



maxima are observed. Fig. 5.4b only displays the motion of the poles for a short time, in order 
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(a) 



(b) 



Figure 5.3: An N = 4, M = 2 Treibich-Verdier solution, with L\ = 4 and L2 = vr. (a) 
The motion of the poles in the fundamental domain. The initial position of the poles 
is indicated by the black dots. The initial time t = was chosen different from the 
collision times t c . The arrows denote the motion of the poles, (b) The KdV solution 

u(x, t) 



not to clutter the picture. The motion of the poles is presumably quasiperiodic in time, as is the 
case for the (4 : l)-solution. It appears that the two poles above (or below) the middle line of the 
fundamental domain share a common trajectory. It is an open problem to establish whether or not 
this is the case. 



5.4 N = 6: two different possibilities. A three-phase solution 

For N = 6, two distinct pole configurations are possible. The first one corresponds to a (5 : 1)- 
solution and results in a two-gap potential of the Schrddinger equation. It essentially behaves 
as the (3 : l)-solution with two more poles added, which also behave as the outer poles of the 
(3 : l)-solution. 

The second configuration is a (3 : 2 : l)-solution, which limits to a three-soliton solution with 
wavenumber ratio (k\ : k% : £3) = (3:2:1). This elliptic solution is a three-phase solution of the 
KdV equation. 

The amplitude ratio of the (3:2: l)-solution is (9 : 4 : 1), which explains why the third phase 



is hard to notice in Fig. 5.£. Animations of the pole dynamics and of the time dependence of the 



(3:2: l)-solution are available at http : //amath-www . Colorado . edu/appm/other/kp/papers| . 
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Figure 5.4: jV = 5: (a) The pole dynamics of a (4 : l)-solution in the fundamental 
domain, (b) The pole dynamics of a (3 : 2)-solution in the fundamental domain, (c) 
The KdV solution u(x,t) corresponding to the pole dynamics in (b). In (a) and (b), 
the initial positions of the poles are indicated. Both solutions are quasiperiodic in 
time. u(x,t) 
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Figure 5.5: The pole dynamics of a (3 : 2 : l)-solution. The black dots mark the initial 
position of the poles; the grey dots mark the position of the poles at t = 0.4. 




Figure 5.6: A (3 : 2 : l)-solution of the KdV equation, corresponding to the pole 
dynamics in Fig. |5.5| 
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